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Abstract. We prove that the independence complexes of some grids have exponential Betti 
numbers. This corresponds to the number of ground states in the hard-core model in statistical 
physics with fermions in the vertices of the grid. 



1. Introduction 

The purpose of this paper is to investigate some topological questions arising from the study 
of supersymmetric lattice models in statistical physics. Suppose one has a finite graph L, which 
in applications is usually a periodic lattice with some boundary conditions. Square, triangular 
or hexagonal grids are the most notable examples. The vertices of the graph can be occupied by 
particles, such as fermions, which satisfy the hard-core restriction: two adjacent vertices cannot 
be occupied simultaneously A configuration of particles which satisfies this assumption is an 
independent set in the graph L. 

Associated with L there is a simplicial complex called the independence complex of L and 
denoted I(L). Its vertices are the vertices of L and its faces are the independent sets in L. It is a 
standard object studied in combinatorial algebraic topology. There is a close connection between 
the simplicial and topological invariants of I(L) and certain characteristics of the corresponding 
lattice model which are of interest to physicists. It is beyond the scope of this paper to discuss 
this relationship in detail; we refer to [10] and we limit ourselves to presenting just the most basic 
dictionary: 

the partition function of L the /-polynomial of I(L), 

the Witten index of L minus the reduced Euler characteristic — 

the number of zero energy ground states the dimension of H*(I(L);Q). 

There has been some very successful work calculating the Witten index, homology groups or indeed 

the complete homotopy type of the independence complex for various lattices, eg. [3J |H1 ES HU 
[131 UM HI]. In this paper we focus on the large-scale picture. Computer simulations of van 

Eerten [3] indicate that for some types of lattices, as their size increases, the number of ground 

states grows exponentially with the number of vertices, that is 

dimH*(I(L)) ~ a v M 

for some constant a depending on the type of the lattice, where v(L) denotes the number of vertices 
in a graph L. This situation is called superfrustration and has interesting physical implications, 
see [5] . Engstrom [5] developed a general method of computing upper bounds for the constant a. 
For the lattices of [3] it gives bounds very close to the values predicted in [3]. 

This paper has two main parts. In the first one we present a method which can be used to 
construct exponentially many linearly independent homology classes in I(L) for graphs L of certain 
type. That proves superfrustration of certain lattices and we give examples based on modifications 
of the triangular lattice. In the second part we prove a generalization of the main result of [5], 
which can sometimes give better upper bounds. 
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Both methods work particularly nicely with one type of lattice studied in [3J [5] : the hexagonal 
dimer, also known as the Kagome lattice (see Figfl]). Under suitable divisibility conditions on the 
height and width we will prove that a graph H of that type satisfies: 

1.02"™ » (2V36)»(H) < dimj^H)) < (14 1 / 36 ■ 2 1 / 6 ) v ^ » 1.21"™. 

We will prove the lower bound in Section [3] and the upper bound in Section [5] The previous upper 
bound of [5] was 2 1 / 3 « 1.26 and the experimental approximation by [3] is 1.25 ± 0.1. 

Our technique for lower bounds produces slightly more than just homology classes: we obtain 
a large wedge of spheres that splits off. For instance, for a suitable lattice H of Kagome type, this 
reads as a homotopy equivalence 

I(M) ~ ( V/ 5 2 "™/ 9 - 1 ) V X 
for some space X. This type of result is proved in Section [I] 

Remark 1.1. Estimations as above can be compared against the absolute upper bound: for any 
graph G we have 

dim H,(I(G)) < (2 2 / 5 )"< G ) « 1.32"^. 
This follows from the results of [T7] ; for another short proof see [T] . 

Acknowledgement. The author thanks Alexander Engstrom, Liza Huijse and Kareljan Schoutens 
for very helpful suggestions. 

1.1. Notation. For a graph G we denote by V(G) the set of vertices and by v(G) its cardinality. 
For any vertex v we write N[v] for the closed neighbourhood of v, that is the set consisting of v 
and all its adjacent vertices. For any set W C V(G) we define A[W] = U ve wN[v]. 

For any simplicial complex A and a subset U C V(K) of the vertices K[U] denotes the induced 
subcomplex of A with vertex set U. The same notation is used for graphs. If H is an induced 
subgraph of G then 1(H) is an induced subcomplex of 1(G). By |A| we denote the number of 
faces in A, including the empty one. The join A * L of two complexes A and L is the simplicial 
complex with faces of the form a U r for a £ K and t £ L. For any two graphs G and H, if G U H 
is their disjoint union, we have 

I(GUH)= 1(G)* 1(H). 

The symbol S k denotes the topological sphere of dimension k. The suspension S K is the join 
K * S°. If e is the graph consisting of a single edge then 1(e) = S°. 

The reduced homology and cohomology groups of A', denoted H*(K), H*(K), are the homology 
groups of the augmented chain, resp. cochain complex of K. Throughout the paper we always 
use rational coefficients and omit them from notation. We have Hi(T, K) = Hj_i(A). There is a 
standard bilinear pairing, denoted (■, ■): 

(•,■) :ir(A)®i^(A)^Q 

given by evaluating cochains on chains. The z-th reduced Betti number is Pi(K) = dim Hi(K) 
and the total Betti number of K is 

p(K) = Y,UK). 

i 

It satisfies /3(£ K) = /3(A). 

Note the empty simplicial complex which has no vertices and a unique face 0. It satisfies 
£0 = S°, so it is good to think of it as S" -1 . It has a single reduced homology group A_i(0) = Q 
and in particular /3(0) = 1. Of course for every non-empty space A we have Hi(K) = for i < 0. 

We do not distinguish between a simplicial complex and its geometric realization. The symbol ~ 
means homotopy equivalence. The reference for other notions of combinatorial algebraic topology 

is pi]. 
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2. Cycles defined by matchings 

A matching of size k is the disjoint union of k edges. An induced matching in a graph G is a 
matching which is an induced subgraph of G. Explicitly, it is a set of k edges of G such that any 
two vertices from distinct edges are non-adjacent in G. 

Definition 2.1. Suppose M is an induced matching in a graph G. A subset a of the vertices of 
M will be called a transversal if it contains exactly one vertex from each edge of M . A transversal 
is dominating if it is a dominating set in G i.e. if N[o~] = V(G). 

Every transversal of an induced matching is an independent set in G. It is dominating if and 
only if it is a maximal independent set. 

If M is a disjoint union of k edges then I(M) = S° * ■ ■ ■ * S° = S k ~ 1 is the boundary of the cross- 
polytope. If M is an induced matching in G then we obtain an embedding S k ~ 1 — I(M) ^ 1(G). 
We will denote by ctM £ Hk-\(I(G)) the image of the fundamental class of S* -1 under this 
embedding (there is a choice of orientations involved, but it does not matter). Explicitly, if 
M = W\), («2 ; ^2)1 ■ • • , (ffc, Wk)} then oim is represented by the cycle 

(1) (M - [vh]) A (H - M) A • • • A ([v k ] - [w k ]). 

If a is an independent set of cardinality k in G then we denote by er v £ C k ~~ 1 (I(G)) the cochain 
which associates ±1 to the two orientations of a and to all other simplices (again, this depends, 
up to sign, on the choice of orientation for a). If a is a maximal independent set (hence a maximal 
face in 1(G)) then this cochain is in fact a cocycle, hence it determines an element of H k ~ 1 (I(G)) 
which we continue to denote tr v . 

Lemma 2.2. If (M, a) is an induced matching with a dominating transversal then q.m an d 17 V 
are nonzero classes in H*(I(G)) and H* (1(G)) , respectively. Moreover 

(cr v ,a M ) = ±1- 

Proof. The last statement holds because c v evaluates to ±1 on exactly one of the simplices in the 
chain representation of au- It immediately implies that both elements are nonzero. □ 

This method of constructing homology classes qm has been used before, also in the context of 
grids [15] where they are called cross-cycles. 

Example 2.3. Consider the cycle Cg with vertex set V(Ce) = {1,2,3,4,5,6} and with two 
matchings with dominating transversals: 

Mi ={(1,2), (4, 5)}, ^ = {1,4}, 

M 2 = {(2,3),(5,6)}, (7 2 = {2, 5}. 

Assuming that the above order defines positively oriented simplices we get: 

(o-y,a Ml ) = 1 for i = 1, 2, (a^,a M2 ) = 0, (o%,a Ml ) = 1. 

It follows that a Mi and au 2 are linearly independent elements of Hi(I(Cq)), hence they are a 
basis, as I(C 6 ) ~ S 1 V S 1 (see [TS]). It means that if M 3 = {(3, 4), (6, 1)} is the third matching of 
this kind in Cq then we must have «m 3 = 010^ + o,2CtM 2 f° r some a\ 1 a^. Applying o~\ and a 2 
to this equation we get —1 = a\ and = a± + ct2 hence ctM 3 = —olmx + olm 2 - 

Remark 2.4. In general, if ai,...,otk £ Hi(X) and 71,..., 7; £ H l (X) then the i-th Betti 
number fii(X) is at least as big as the rank of the I x k matrix with entries (7^, a t ) for 1 < s < Z, 
1 < t < k. 

3. Hexagonal dimer and related grids 

The hexagonal dimer or the Kagome lattice is the lattice obtained from the triangular lattice 
by erasing every other line in each direction in the way shown in Figjl] It is invariant under the 
translations by (2,0) and (0, V&). Let i„ ]m denote the quotient of that lattice by the action of 
the translation group generated by the vectors n ■ (2, 0) and m ■ (0, v3)- It has v(H rajm ) = 3nm 
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Figure 1. The hexagonal dimer lattice, its tiling with large hexagons and two 
types of tiles. 

vertices and if one additionally assumes that 6\n and 4|m then it can be tiled with large hexagons 
in the way shown in Fig[l] 

In fact there is no harm forgetting about n and m. Let EI be any quotient of the hexagonal 
dimer lattice (i.e. its finite portion with cyclic boundary conditions) with the property that it 
can be covered in this way by the large hexagons. Then the number of those hexagons is always 
u(H)/36 and we have the next result. 



Proposition 3.1. For a hexagonal dimer lattice EI which admits a tiling as in Fig^ we have 

/3(/(H)) > (2 1 /36)^(H)^ 

Proof. Let k = w(H)/36 be the number of large hexagons in the tiling of EL We are going to 
construct 2 k linearly independent elements in _f/8fc-i(-f(EI)). That proves the claim since then 
/?(7(H)) > 2 k = 2"( H )/ 36 . 

Consider the tiles A and B of FigJT] where they are shown with an induced matching (thick 
edges) and its dominating transversal (thick vertices). Note that for any placement of A and B 
in place of the large hexagons in the grid we obtain a valid induced matching with a dominating 
transversal in EI. Indeed, there are no edges directly between the tiles, so the matching is induced 
and each vertex located between the tiles is adjacent to one of the transversal vertices on the outer 
cycle of a tile. We therefore have 2 k homology (and cohomology) classes in /(H). 

To prove linear independence we need some notation. Suppose that the big hexagons in the 
grid are labeled 1, . . . , k in some order. For any sequence s — (si, . . . , Sk) of letters A, B let M(s) 
and cr(s) denote the matching and transversal obtained in the above construction by placing the 
tile of type Si in the i-th spot for i = 1, . . . , k. Then for any two sequences s and t we have 

, , . v \ f if ti = B and Si — A for some i 

(<r(t) ,a M{s) ) = | ±1 otherwise 

Consider a 2 fe x 2 fc matrix with rows and columns indexed by sequences in {^4, B} k in lexicograph- 
ical order, where the entry in column t and row s is (o~(t) v , aM( s ))- By the last observation that 
matrix has ±1 on the diagonal and above the diagonal (where t >lex s), so it is of full rank. 
Therefore all cim( s ) are linearly independent. □ 
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Figure 2. Two other lattices derived from the triangular lattice and their sets of tiles. 



Similar results can be obtained for some other lattices derived from the triangular lattice. 

Proposition 3.2. For d = 3,4 let A^ be any finite quotient of the lattice obtained by removing 
every d-th line in each direction from the triangular lattice, such that the new lattice admits a 
tiling as in Fig^ Then we have 

/3(I(A 3 )) > (2 1 /8)«(^3) ; /3(/(A 4 )) > (2 1 /45)^(A 4 )^ 

Proof. The proof is similar: A3 has k = v(A.s)/8 tiles and A4 has k = t>(A4)/45 of them. In 
each case we construct 2 fc induced matchings with dominating transversals using the tiles A, B 
of appropriate type (Fig[5]) and prove linear independence of the resulting homology classes as 
before. □ 



4. HOMOTOPICAL SPLITTINGS 

We will now show for completeness that induced matchings with dominating transversals in G 
correspond to sphere wedge summands in the homotopy type of 1(G). This is stronger than just 
saying that they define a nonzero homology class. The results of this section are not needed in the 
previous calculations, where analyzing the cohomology-homology pairing is sufficient (and easier). 

Lemma 4.1. Let (M,a) be an induced matching of size k with a dominating transversal in G. 
Denote by j : I(M) <^-> 1(G) the embedding induced by the inclusion i : M <^-> G. Then there is a 
homotopy equivalence 

1(G) ~ I(M) V C(j) 
where C(j) is the homotopy cofibre of j . 

Proof. Since a is a maximal face of 1(G) we can remove it and form the simplicial complex I(G)\o~ 
whose geometric realization is obtained by removing the interior of a. Then we have a cofibre 
sequence 

da 1(G) 1(G) -> S (da) -> • • • 

i.e. 1(G) is the homotopy cofibre of the inclusion da (1(G) \ a). But this inclusion is null- 
homotopic since it factors through the space I(M) \ a, homeomorphic to the disk D k ^ 1 . As a 
consequence 

1(G) ~ (1(G) \ a) V S (da) = (1(G) \ a) = (1(G) \ a) V S^ 1 . 
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Now the fact that I(M) \ a is a contractible subcomplex of 1(G) \ a implies that 

1(G) \ a ~ (1(G) \ a)/(I(M) \ a) = I(G)/I(M) ~ C(j). 
Together with S k ~ 1 — I(M) this ends the proof. 



□ 



To analyze our situation further we need to understand homotopy cofibrcs of maps between 
independence complexes. 

Lemma 4.2. Suppose i : H '—t G is an inclusion of an induced subgraph. Let C(i) denote the 
graph obtained by adding to G a new vertex adjacent to all the vertices ofV(G) \ V(H). Then 
I(C(i)) is the homotopy cofibre of the induced inclusion 1(H) 1(G), i.e. there is a cofibre 
sequence 

1(H) <—> 1(G) <^> I(C(i)) S 1(H) • • • . 

Proof. By definition the homotopy cofibre of 1(H) ^ 1(G) is obtained by attaching to 1(G) a 
cone over the subspace 1(H). The space I(C(i)) is obtained in precisely the same way. □ 

Now we can prove the main result of this section. 

Proposition 4.3. Suppose (Mi,o~i), for i = l,...,k, is a sequence of induced matchings with 
dominating transversals in a graph G, such that for every i < j we have 

(*) *j \ V(M i ) 0. 

Then there is a homotopy equivalence 

k 

/(G) ~ (V/ I (Mi)) VI 
i=i 

for some space X . 



Proof. Let i\ : M\ H G be the inclusion of the first matching. Then, by Lemmas 4.1 and 4.2 we 
have a splitting 

1(G) ~ /(Mi) V I(C(h)). 

If k = 1 then we are done. Otherwise note that the pairs (Mj , o-j ) for 2 < j < k define a sequence 
of induced matchings with dominating transversals in the graph C(i\). Indeed, to build C(i\) we 
did not add any edges within G itself, so the matchings are still induced. Moreover, by Q every 
set a-,- contains a vertex which is not in V(Mi), hence it is adjacent to the new vertex of C(i\), 
which means o~j is a dominating set in C(i\). The condition ([*]) still holds, so by induction 

k 

I(C(h))~ \/ I(Mj)VX 

3=2 

and that completes the proof. □ 

Remark 4.4. The last result applies to all the situations of the previous section, once one orders 
the matchings in the lexicographical order of their defining {A, B}-words. In particular we get a 
splitting of /(H) which includes a wedge sum of exponentially many spheres. 

5. Upper bounds 

In this section we improve, for the hexagonal dimer grids H, the upper bound of [5]. We 
develop a more general result, which is modeled entirely on the technique of [5] with just two 
small improvements. Firstly, it avoids discrete Morse theory and relies just on the Betti numbers 
and homotopy. Secondly, it allows arbitrary graphs in a place where [5] requires a forest (see 



Cor 5.41 



We first need some extra notation. Suppose If is a simplicial complex with a fixed splitting of 
the vertex set into two disjoint subsets U and W (V(K) = U U W). For every simplex a € K[U] 
we denote by lk w a the subcomplex of if [VF] consisting of those simplices r G ZT[W] for which 
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r U a G K. By st^/a we denote the subcomplex of K consisting of those r G if for which r n £/ C a 
and r n W G lkv^c. Clearly st^tr = (lk^o - ) * c 

Example 5.1. If a = then we always have lk^0 = stw0 = 2<r[W]. If U = {v} then the 
complexes Ykw v and stv^f coincide with the usual link and star of v in K. The star stiver is 
always contractible when a ^ 0. 

Lemma 5.2. Suppose K is a simplicial complex with a vertex partition UU W as above and such 
that for every a G K[U] we have /3(lkn/<j) < B. Then 

P(K)<B-\K[U}\. 

Proof. Denote D = \K[U]\. Fix any ordering = cr > ■ ■ ■ > a D-\ of the simplices in -K"[f7] such 
that every simplex is preceded by all its faces. For < I < D define subcomplexes 

l 

Fl = U st W°i- 

i=0 

Then JC[W] = Fq C Fi C • • • C Fd_i = K is an increasing, exhaustive filtration of K with 
quotients: 

Since each of those quotients has total Betti number at most B, the E\ page of the homology 
spectral sequence associated with this filtration has dimension at most B D. The spectral sequence 
converges to H*(K), hence the result. □ 

Lemma 5.3. Suppose U C V(G) is a vertex set with the property: For every simplex a £ I(G[U]) 
the total Betti number of 

I(G\(UUN[a])) 

is at most B. Then 

P{I{G)) <B-\I(G[U})\ < B-2W. 

Proof. Consider K = 1(G) with vertex set partitioned into U and W — V(G) \ U. Then for every 
simplex a G I(G[U}) we have precisely 

lk w a = I(G\(U{JN[a])), 

so Lemma |5.2| applies. □ 



Corollary 5.4 ([S]). If U C V(G) is a vertex set such that G\U is a forest then f3(I(G)) < 
\I(G[U])\. 

Proof. If G \ U is a forest then so is G \ (U U iV[cr]) for every a. Since by [1] the independence 
complex of a forest is either contractible or homotopy equivalent to a sphere (possibly S 1-1 ), we 
can apply Lemma |5.3| with B = 1. □ 



Fact 5.5. For any iwo topological spaces X and Y 

Proof. This follows from the formula for the reduced homology of the join (eg. [20, Lemma 2.1]), 
which for rational coefficients reduces to 

H k (X*Y)= Hi(X) ® Hj(Y), k>-l. 



i,j>-i 

i+j=k-l 



□ 



Proposition 5.6. IfM. is any hexagonal dimer lattice which can be tiled as in Fig\l\ th 

P(U) < (14V36. 2 V6)HH)_ 
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Figure 3. The induced subgraph G of the 30-vertex tile with /3(I(G)) = 14. 

Proof. Recall that H contains k = w(H)/36 large hexagons. Let U consist of those vertices of H 
which are not covered by the tiles. Then \U\ — 6k — u(H)/6 and the graph H \ U is a disjoint 
union of k hexagonal tiles with 30 vertices each. The graph M.[U] has no edges. 

Now suppose that a is any subset of U. Then the graph H \ (U U N[a]) is a disjoint union 
of k graphs, each of which is an induced subgraph of the 30-vertex tile. More precisely, it is 
an induced subgraph obtained by removing the neighbourhood of some subset of the vertices of 
U which surround that tile. There are 2 12 = 4096 graphs that can arise in this way, with only 
217 isomorphism classes |21| . and the homology of their independence complexes can be easily 
calculated by a computer [TS] . It turns out that for each of those graphs the total Betti number 
is at most 14 (the graph which attains maximum is shown in Fig{5]). Using Fact 5.5 we get that 
the total Betti number of H \ (U U N[a]) is at most B = 14 fc . 

Lemma |5.3| now gives the conclusion: 

/3(/(H)) < 14 fc • 2 |c/| = 14"( H )/ 36 • 2 v{w)/& . 

□ 

Remark 5.7. If the dimensions of the grid H do not allow it to be tightly tiled with the large 
hexagons then Prop |5.6| still holds asymptotically. This is because one can pack H with hexagons 
leaving just a region of size proportional to the perimeter of the grid. The wasted vertices can 
then be added to the set U . 

Remark 5.8. Using exactly the same technique one proves the following counterparts to the lower 
bounds of 13.21 

/3(/(A 3 )) < (2 3 / 8 f (A3) , /3(I(A 4 )) < (10 1/45 • 2 1 / 5 ) v{ ^\ 

6. Concluding remarks 

It would be most interesting to prove similar lower bounds for the classical triangular and 
hexagonal lattice, for which superfrustration is also predicted [3i. In those cases, however, the 
present methods do not seem to work, and it is not clear if exponentially many homology classes 
should be given by embedded spheres or if more complicated constructions are necessary. Recent 
results in this direction include [TT], where "long and thin" triangular lattices of size c x n for 
small constants 2 < c < 7 are investigated. Also, the constructions of Jonsson [T5] can be adapted 
to show that for the hexagonal grids of suitable sizes c x n with fixed c the number of ground 
states is exponential in n. 

More generally, it would also be interesting to know what aspect of regularity is responsible for 
superfrustration in arbitrary lattices. 
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